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ABSTRACT 


Coupled  heterogeneous  electron  transfer/homogeneous  reaction  of  the  type 


A  +  e"  ♦  B  Ej* 

0  +  e"  ♦  R  E2* 

mA  +  R  ♦  AmR  K 

occur  in  the  electrooxidation  of  certain  anilines,  and  in  association 
reactions  of  hydrocarbon  anions  with  metal  cations.  If  the  adduct  A_R  is 
electroinactive,  the  resulting  vol tammogram  is  distorted  in  a  way  which  is 
characteristic  of  the  value  of  m.  The  derivation  of  the  theoretical  fluxes 
for  these  reactions  and  calculated  normal  potential  pulse  vol tammograms  are 
presented.  The  results  are  derived  from  both  limiting  flux  and  digital 
simulation  considerations. 


INTRODUCTION 


We  consider  in  this  report  a  coupled  sequence  of  heterogeneous  electron 
transfer  and  homogeneous  chemical  reactions  in  which  the  chemical  reactions 
involve  the  species  of  two  different  redox  pairs. 


A  +  e”  *  B  Ef 


Cl] 


0  +  e"  x  R  E2*  ’  [2] 

kf  [A]*m  kf 

mA  +  R  rirt  a_r  k'  *  — r - 1  [3] 

kb  b 

where  [A]*  is  the  bulk  concentration  of  A. 

Reactions  second  order  and  higher  in  A  and/or  R  are  also  considered. 
Examples  of  these  reactions  Include  the  reaction  of  hexae thy 1  benzene  cation 
radical  with  hexamethyl benzene  (1)  and  other  hydrocarbons  (2),  and  the 
compounds  in  part  II  (3)  of  this  study,  for  example: 


Li+  +  e"  j  Li(Hg)  Ej*  [4] 
Anthracene”  +  e"  ♦  Anthracene*  E2*  [5] 
3Li+  +  Anthracene*  ♦  Li3(Anthracene)+  [6] 


and 


Na+  ♦  e‘  ♦  Na(Hg) 


[7] 


Fluoranthene  +  e"  ♦  Fluoranthene" 


[8] 


E2 


Na+  +  Fluoranthene'  ♦  Na (Fluoranthene) 


C9] 


The  voltammetry  of  these  systems  exhibits  rather  spectacular  properties  such 
as  depression  of  faradaic  currents  to  zero  and  in  some  circumstances  formation 
of  “pseudo"  peaks  and  dips.  These  effects  in  the  cyclic  voltammetry  of  the 
hexaethyl benzene  oxidations  are  explained  by  reaction  with  other  electroactive 
hydrocarbons  present  wherein  the  second  hydrocarbon  concentration  is  therefore 
depleted.  As  a  result,  the  flux  of  A  to  the  electrode  is  reduced,  resulting 
in  current  depression. 

Unfortunately,  analytical  solutions  for  the  current  and  concentration 
profiles  are  precluded  by  the  non-linear  form  of  the  diffusion-kinetic 
equations  which  describe  the  reaction  sequence  [l]-[3].  Some  predictions  of 
the  general  behavior  can  be  derived,  however,  and  digital  simulation  of  the 
system  is  quite  straightforward,  and  provides  details  of  the  expected 
responses  as  a  function  of  the  usual  parameters. 


Theon 


It  Is  assumed  that  (a)  B,  R,  and  A^R  are  electroinactive  in  the  potential 
region  considered  and  (b)  AmR  is  chemically  stable.  The  experimental 


technique  considered  will  be  normal  pulse  voltammetry.  This  technique  is 


useful  in  the  study  of  systems  where  species  A  or  0  is  an  alkali  metal  ion 


(part  II  (3)).  The  alkali  metal  is  then  soluble  in  and  deactivated  chemically 


by  the  mercury  drop,  and  is  stable  against  interfering  reaction  with  other 


system  components.  In  the  experiment,  increasingly  larger  potential  pulses 


are  applied  to  new  large  mercury  drop  electrodes.  The  diffusion  to  such 


electrodes  may  be  approximated  as  linear.  The  current  transient  Is  sampled  at 


identical  time  intervals  on  each  drop.  The  time  interval  is  usually  chosen  so 


that  charging  current  magnitudes  have  decayed  to  values  well  below  the 


faradaic  component  so  that  their  effect  on  the  total  current  is  minimized. 


The  simulation  of  such  an  experiment  Is  thus  developed  by  considering  a 


succession  of  chronoamperometric  responses,  each  to  slightly  higher  potential 


pulse  values  at  a  planar  electrode.  Corrections  to  the  responses  for 


sphericity  effects  for  the  cases  considered  herein  proved  to  be  less  than  3S 


in  the  worst  case,  and  will  be  neglected. 


The  partial  differential  equations  describing  the  system  [7]-[9]  are 


^  *  da  ^  "  mkfCA3mCR]  +  mkb  [AmR] 


3[0]  ,  Dn  *  [0] 

t r  Do - 


4r  *  °r  “T  “  m  kf[A:j,nCR:i  +  m  “bty^ 


3  EAR]  a^ramol  _ 

sr-  *  V  4 


where  T,  X,  and  are  time,  distance  from  the  electrode  surface,  and 
diffusion  coefficients  of  the  i  th  species,  respectively. 

The  conditions  for  the  experiment  are  summarized  as: 


Initial  (T  *  0): 


All  X:  [A]  =  [A]*,  [0]  =  [Of,  [R]  «  [A^R]  =  0 

[14] 

¥ 

Boundary: 

T  >  0,  X  *  0:  [A]0  *  [B]°  exp  [F(E— E1° )/RT]  * 

[15] 

[0]°  *  [R]°  exp  [F(E-E2°)/RT]  «  62 

[16] 

(o  represents  surface,  X  *  0) 

D  r  -n  f  ^  C  B  3  \ 

DA(“TT}X  =  0  VTT'x  *  0 

[17] 

D  1 )  «  -n  (3[R]v 

do'tt"jx  *  o  V“T r]x  «  o 

[18] 

1  "  FA  Cda(^}x  *  0  +  Do(^}x  -  0] 

[19] 

T  >  0,  X  «  »: 

idem  [14]  [20] 

The  Cottrell  currents  at  any  potential  for  the  processes  [7]  and  [8]  in 
the  absence  of  any  homogeneous  reactions  are  given  by 


where  A  is  the  electrode  area. 


and 


10  =  FAD,,1'2  [03*1(1  ♦  (^)1/2i2).1/2T1/2]-1 

R 


[22] 


The  sum  of  these  two  currents  represents  the  total  Cottrell  current  at  any 
potential.  The  actual  current  for  the  kinetic  system  will,  however,  be 
affected  by  the  decrease  In  [A]  near  the  electrode  due  to  the  chemical 
reaction  [9]. 

We  define  the  following  dimensionless  parameters: 


t  *  jr,  r  *  The  chronoamperometric  pulse  width  in  seconds 


r*  C0  „  [0]  „  [A] 

[A]  ■  CA,  a  =  ,  Cq  *  ,  CA  ■ 


Cd  s  -t-  ,  Ca  d  s  — x 

C.  Am  CA 
A 


x  = 


<DAr> 


m 


kfCA 


Bfm 


e  8  kfCA  +  *b  y  ‘TTT 


[23] 


[24] 


[25] 


[26] 


Equation  [10]  becomes: 


3^a  ^A^A  3 


=  D.r  ^f<cJ)m+1CAmCR  +  mJcbCAC* 


JA  3X2 


*  7^  -  ”¥<<>’ CR  ♦  "Vc*  « 

a  a  m 

3^C 

‘  1*  '  CI  CR  *  V3 


Similarly, 


3C0 

•% 

TT  * 

Q - 7 

3X^ 

3CR 

r3% 

TT  = 

3  X2 

aDAmR 

m 

32C 

TT 


*  s 


-  t[k(c;i"-1  c;  cR  *  cVR] 


*  ’ftK(cI)m'1  c;  cR  *  cA_R: 


m 


The  new  initial  and  boundary  conditions  become 


Initial  (t  *  0)  * 

AH  X:  CA  -  1,  Cp  -  4  .  <*  ■  C^r  •  0 


Boundary: 
t  >  0,  x  *  0:  C? 


C°  exp  [n,  F(E-E,°)/RT] 


[36] 


Cj  *  Cj  exp[n2F(£-E2°)/RT]  *  62 

l7TT)x  =  o 5  V  "u(rr}x  *  o 

3C  3Q 

q(TJT)X  =  0  =  *0  s  ~r[Trh  =  0 
t  >  0,  X  •  :  idem  [34] 


[37] 

[38] 

[39] 


The  total  electron  flux,  in  dimensionless  units,  follows  from  [19]: 

♦  -  *A  +  *0  "  1(J-)1/2(FAC*)‘1  [40] 

UA 

The  dimensionless  Cottrell  fluxes  for  the  two  electron  transfer  processes 
are  given  by  (see  [21]  and  [22]) 


A  (Cottrell 


j-[(l*  u'1/2  SjMwt)172]"1 


*0  (Cottrell )  *  C(1  +  <?)l/2 


[41] 

[42] 


for  any  potential  pulse  magnitude.  When  the  processes  [7]  and  [8]  are 
diffusion  limited,  [41]  and  [42]  reduce  to 


*  -  c*l  M/2 

*A  (Cottrell,  D.L.)  " 

hi  -  a  rW/2 

*0  (Cottrell,  D.L.) 


[43] 

[44] 


so  that  the  Cottrell  currents  are 


F<  ( 


;A,l/2 


itT 


[45] 


1-o 


FA  aC 


*,qVl/2 

AMTtF 


[46] 


Equations  [30]  -  [33],  taking  the  conditions  [34]  -  [39]  into  account 
were  solved  by  finite  difference  and  orthogonal  collocation  techniques  to 
provide  the  concentration  profiles  necessary  to  deduce  the  current  response 
[40]. 

A.  value  of  m  in  equation  [9]  may  be  determined  analytically  when  (a)  both 
electron  transfers  are  diffusion  limited  (high  values  of  -E)  and  (b)  kf  is 
very  large.  The  boundary  conditions  then  stipulate  that  C®  and  Cj  =  0.  Two 
limiting  cases  will  then  exist  that  describe  the  total  electron  flux  *e  to  the 
electrode. 


<1>  Assume  nn|>Q  > 


[47] 


If  k  -*•  »  then  CA  will  be  depressed  to  zero  concentration  wherever  mCR  > 
C A.  Thus  the  value  of 


3Ca 

(“57)X  =  0  [48] 


must  be  zero,  since  C®  is  also  at  zero  (Figure  1).  Hence  A  will  not 
contribute  to  the  electron  flux  and 


[50] 


i  =  FA*0  •  FA^fl^)1'2 


<2>  Assume  m*0  <  ^ 


If  k  •,  (Figure  2)  the  flux  of  A  is  reduced  by  m  times  the  flux  of  0, 


♦e  ‘  *0  +  <**  ’  V  ’  *A 


i  ■  FA*,  =  F»C;  (^)1/2 


If  the  electron  fluxes  given  by  [49]  and  [52]  are  normalized  to  i»0  ,  we 


*e  *0 

<1>  For  J-  -  ji  -  1 

0  *  wo  wo 


<2>  For  *,>*„,  ^  ^ 


It  will  be  of  interest  later  to  consider  cases  where  the  number  of 


electrons  transferred  to  the  0  species  is  greater  than  1,  i.e. 


t)2  >  1 


msssmmimsssBsm 


Therefore,  proceeding  as  above,  the  general  case  for 


I 

>v 


K- 


r-.v. 


V\ 

& 


IS 


A  +  e"  *  B 


0  *  n2e“  *  R 


mA  +  R  ♦  AmR 


[57] 

[58] 

[59] 


leads  to 


<1> 


s  n*  for  r4u  >  ♦. 


[60] 


and 


<2> 


♦e  n2^o  +  ^A  “  "V  _  (n2  '  *0  +  ^A 

’  ^  C 


[61] 


C*  D 

1  LA  UA 

(n2  -  m)  + - jty  B  *n2  "  m)  +  t — n  177 


[62] 


oq 


C0  D0 


for  >  m^Q 


[63] 


ve  1 

Hence  a  plot  of  the  normalized  current  -r-  vs.  - yT7 

*0  a  q 


has  the  form  shown 


in  Figure  3.  The  point  of  intersection  of  the  two  linear  portions  will  be 
when 


^0  =  ♦a 


[64] 


At  that  point 


The  y  intercept  for  the  extrapolated  -4-  line  for  case  <2>  will  be  at 


(m  -  n2) 


[66] 


It  should  be  pointed  out  that  there  are  other  general  shapes  for  the  plot  In 
Figure  3,  depending  on  the  values  of  n2  and  m.  The  rules  for  the  analysis 
treatment  are  still  the  same. 


SIMULATION 

(A)  FIRST  ORDER  KINETICS  in  CA  and  CR  (i)  E2*  >  Ej’ 


From  equation  [26],  we  have  the  relation  between  the  real  rate  constants 
equilibrium  constant,  dimensionless  parameters,  and  pulse  length: 


yK 


C?  rm 
A 


y 


cJrn*f 


T 


[67] 


*  Ti; 


T  mkk 


[68] 


The  results  presented  in  this  section  will  be  for  a  measurement  made  at 
10S  of  the  total  chronoamperometric  pulse  time.  Figure  4  shows  the  simulated 
normal  pulse  voltammograms  for  a  fixed  concentration  of  species  A  as  a 
function  of  the  0  concentration.  In  this  case  y  *  K  *  0,  so 


that  and  are  given  by  the  Cottrell  relations  [41]  and  [42].  The 


diffusion  currents  represent  simple  uncoupled  charge  transfers.  The 
concentration  increments  of  the  0  species  are  equal  at  0.20  C Ab  units,  and 
range  (bottom  to  top)  from  0.0  to  4.0  CAb  units.  We  assume  DA  *  D0  for 
simplicity. 

Figure  5  represents  the  same  system  now  perturbed  by  the  kinetic  reaction 

A  ♦  R  AR  K  *  [69] 

kb  b 


The  reaction  [69]  is  first  order  in  each  reactant  A  and  R.  In  the  Figure,  y  * 
4.0  and  K  *  1.0.  It  is  seen  that  at  all  concentrations  C0*  >  1.0  CA*,  the 
incremental  change  in  the  flux  at  high  negative  potentials  Is  precisely  the 
same  as  in  Figure  4.  The  absolute  value  of  the  total  flux  is  decreased  by  an 
amount  equal  to  *A.  At  all  concentrations  Cq*  <  1.0  CA*  ,  it  Is  seen  that  the 
Incremental  changes  In  total  flux  at  high  negative  potentials  have  been 
reduced  to  very  small  levels,  but  are  still  equal.  The  absolute  values  of  the 
fluxes  lie  between  *  ^A  +  0*2  4»A  and  ♦  *  *A.  In  addition,  peaks  and  dips 
have  appeared  in  the  voltammetry  in  the  region  of  the  second  wave  and 

plateau.  The  height  and  depths  of  these  features,  as  well  as  the  potentials 

of  the  critical  points  are  functions  of  y,  K  and  the  concentration  ratio. 

As  the  values  of  y  and  K  are  Increased,  it  Is  found  that  the  peaks  and 

dips  level  out,  and  the  fluxes  for  the  CQ*  <  1.0  CA*  decrease  to  a  limiting 

value  of  *  4»a,  as  is  called  for  by  equation  [52].  For  CQ*  >  1.0  CA*  ,  It 

Is  found  that  ♦  *  ♦q  as  is  called  for  by  equation  [49].  It  is  pointed  out 

*  * 

here  that  for  C0  >  1.0  CA  ,  the  results  for  the  limiting  at  high  negative 

potentials  are  the  same  for  all  values  of  y  >  2.0  and  K  >  1.0. 

It  is  apparent  from  equations  [49]  and  [52],  for  the  first  order  case 


than  that  corresponding  to  4>A  (Cottrell). 

(ii) .  E2*  -  Ef 

The  results  here  are  virtually  the  same  as  for  the  previous  case,  except 
that  there  is  only  one  voltammetric  wave  (Figure  6).  The  absolute  value  for 
the  electron  flux  for  large  values  of  y  and  K  are  again  given  by  [49]  and 
[52]. 

(iii) .  E2*  <  Ef 

In  this  case  (Figure  7),  0  is  reduced  before  A.  Therefore,  it  is  the 
products  of  the  first  electron  transfer  that  are  involved  in  the  kinetic 
reaction  instead  of  the  reactant.  The  diffusion  controlled  portion  of  the 
first  voltammetric  wave  is  always  only  determined  by  the  flux  of  0  to  the 
surface  and  will  therefore  now  be  subject  to  depression  by  the  kinetic 
reaction.  The  height  of  the  second  electron  transfer  wave  will  be  strictly 
determined  by  the  homogeneous  kinetics,  and  again  Is  given  quantitatively  by 
[49]  and  [52]  for  polarograms  for  this  case.  The  half  wave  potentials  are 
shifted  by  amounts  governed  by  the  ratios  of  the  surface  concentrations  of 
each  redox  pair,  just  as  they  are  In  the  other  cases. 

(B).  Multi-Order  Kinetics  in  A  and  R 

The  same  general  conclusions  can  be  made  for  this  case  as  were  made  for 
(A)  1,  11,  and  111  above.  For  high  values  of  y  and  K,  the  diffusion  limited 
current  at  high  negative  potentials  is  governed  by  the  relations 


(see  [60]  and  [61]). 

An  Important  feature  here  Is,  that  when  and  m  >  n2  ,  It  Is 

possible  ( cf  m  •  3,  n2  *  1,  s  0.31,  i|/q  =  0.1  )  for  the  net  current  flux 
to  be  less  than  <|>A.  Hence  plateaus  at  levels  less  than  that  of  ^>A  are 
observed  (Figure  8). 

Results  for  the  reaction 

3A  +  R  ♦  AjR 

with  CAb  =  C0b  as  a  function  of  CAb  is  shown  in  the  Figure  to  demonstrate  the 
diversity  in  possible  voltammetric  shapes  for  the  same  system.  The  kinetic 
parameters  here  are  y  *  4.0,  K  ■  1.0  and  Ej*  <  E2*.  At  high  concentrations, 
two  maxima  and  a  minimum  are  produced  before  diffusion  limiting  flux  Is 
attained.  At  lower  concentrations,  one  maximum  is  seen,  and  at  even  lower 
concentrations,  there  are  no  inflections.  An  example  of  this  type  of  behavior 
is  presented  in  part  II  of  this  series  (3). 
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FIGURE  LEGENDS 

Figure  1.  Typical  concentration  profiles  for  the  case  rt?0  >  *A,  k  *  •> 


Figure  2.  Typical  concentration  profiles  for  the  case  m*0  <  *A,  k  *  » 

% 

Figure  3.  Plots  of  the  normalized  current  function  vs.  the .normalized 
bulk  concentration  -  diffusion  coefficlentOratio  (a  q)"A/z  for: 
(A)  m^Q  >  fA;  and  (B)  rtf  <  ?A. 


Figure  4.  Simulated  normal  pulse  voltammograms  for  two  electroactive 
species  A  (first  wave)  and  0  (second  wave)  with  no  kinetic 
complications.  Bottom  vol tammogram  Cq  =  0.0.  In  each 
successive  vol tammogram,  the  concentration  of  0  is  incremented 
by  0.2  CA  units.  The  top  curve  corresponds  to  Cq  *  4.0  CA.  The 
diffusion  coefficients  are  assumed  equal. 


Figure  5.  As  Figure  4  with  y  *  4.0  and  K  =  1.0. 


Figure  6.  Simulated  normal  pulse  voltammograms  for  the  case  Ej *«E2* .  The 
concentration  of  CA  Is  held  fixed  and  the  concentration  of  C0  is 
varied  (bottom  to  top)  from  1.0  to  3.8  C*  in  increments  of  0.4 
CA.  The  superimposed  curve  (for  clarity)  is  for  no  kinetic 
reaction  in  a  mixture  of  A  (CA  *  1.0)  and  0  (Cq  *  4.0  CA). 


Figure  7.  Simulated  normal  pulse  voltammograms  for  the  case  E2*  <  Ej*. 

The  concentration  of  CA  is  held  fixed  and  the  concentration  of 
Cq  is  varied  (bottom  to  top)  from  1.0  to  3.8  CA  In  Increments  of 
0.4  CA.  The  superimposed  curve  for  clarity,  Is  for  no  kinetic 
reaction  in  a  mixture  of  A  (CA  *  1.0)  and  0  (Cq  *  4.0  CA). 


Figure  8.  Simulated  normal  pulse  voltammograms  for  the  case  of  multi - 

ordered  kinetics;  3A  +  R  ♦  A3R,  CA  *  Cq,  y  »  4.0,  K  *  1.0,  Ej*  < 
E2*.  From  bottom  to  top:  concentrations  *  0.25,  0.50,  0.75, 
1.0,  1.5,  2.0,  2.25,  and  2.5. 
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